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We showed that various expressions of stress in different 
models of fluids and in different forms of applying Newton's 
Second Law can be chosen, however there is always one New- 
ton's Second Law which restricts and provides flexibility in 
these expressions. When fluids are regarded as being made 
of cells from a microscopic point of view, temperature force 
(gradient of internal kinetic energy) can also drive them flow. 



While the pressure is the normal force acting on a sur- 
face per unit area, the stress means the total force per 
unit area essentially. However, for a given system of par- 
ticle i from 1 to n with volume £1, the stress has been 
proven or widely used as 



7T = 7Tfc + TTf , 

where the kinetic energy term is 
1 - 

8=1 

and the regular force term is 
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with mass rrij and position vector r$ for particle i, and 
force fj—>i acting on particle i by particle j. As the kinetic 
energy term seems not representing any forces directly, 
Zhou recently studied previous derivations for it and con- 
cluded that the kinetic energy term was not physically 
grounded and should be deleted from the stress expres- 
sion [1]. In our recent work [2], we found that the kinetic 
energy term can be explained by collision forces between 
particles and walls surrounding the system in a statisti- 
cal way. As there has been no doubt about the regular 
force term in principle and we have discussed it in detail 
in our previous work for systems with periodic boundary 
conditions [2], we will mainly discuss the kinetic energy 
term in views of macroscopic and microscopic motions 
based on Newton's Second Law below. 
Let us repeat Newton's Second Law first, 
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Fi, (i = 1,2, - • -, n) 



(4) 



where Fj is the net force acting on particle i. Even 
for this extremely well known law, we still want to 



emphasis that it means the rate of the "regular-force- 
induced" change of the moment is equal to the regular 
force essentially. For one particle, this is never a prob- 
lem. For many particles, if we simply add Eq. (4) for all 
particles and get 



(5) 



no problem either. Now let us consider a very short 
period of time At = t% — fi > (from t\ to t% ) for 
a system containing the above n particles but defined 
with a fixed space boundary. Suppose an extra parti- 
cle, particle i + 1, enters the system during this period 
even with an unchanged velocity r n +i, and no particles 
move out of the system. If we calculate the momentum 
of the system at t\ and t% as Pi = X^j=i m i^i i^i) an d 
P2 = m n+1 r n+1 + YZ=i ( i2 ) respectively, and get 
the average rate of the momentum increase of the sys- 
tem for this period as 



At 



Pi r„ + i 



Er 4 (t 2 ) - 1 
mi : — 
At 

i=l 



(6) 



we need an additional imagined force f = m n _|_ir n+ i/At 
acting on the system in order to have Newton's Second 
Law satisfied as 
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(7) 



which is equivalent to Eq. (5). The force f can be imag- 
ined as the interaction between the spaces separated by 
the boundary area penetrated by particle n + 1, as dis- 
cussed in the TRANSPORT OF MOMENTUM section 
in our previous work [2]. So we have two choices. One is 
to keep that the rate of moment change in Newton's Sec- 
ond Law is the rate of the "regular- force-induced" change 
of the moment. Then we have Eq. (5) as Newton's Sec- 
ond Law with the imagined force f never needed. The 
other is that the momentum change due to any reasons, 
including particles' passing system boundary, is accumu- 
lated, then we have Eq. (7) as Newton's Second Law, 

and the imagined forces f must be considered. These 
physical views for particle systems are also applicable in 
complicated cases of continuous medium/fluid, for which 
forms of stress are usually derived and extended to more 
general cases. 
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In 1950, Irving and Kirkwood derived their formula 
for stress similar to the internal stress below by apply- 
ing statistics onto hydrodynamics [3]. We will apply 
Newton's Second Law in hydrodynamics and limit statis- 
tics to necessary terms only, and will deal with various 
stress definitions. It will be convenient for us to start 
our derivation by repeating some of their work. Let us 
consider a continuous fluid with mass density p (r; t) and 
"local velocity" u (r; t) at the point r and time t. We now 
imagine w to be a space region somewhere in the interior 
of the fluid with fixed boundary S. Since w is fixed, the 
rate of change of momentum in it is 

Rc = | J p(r;t)u (r; t)dw = J -| [pu] dw. (8) 

The rate of flow of momentum out through the surface 
of lo is 

Rf — (f> p(r;t)u (r; t) u (r; t) ■ dS — / V r • [puu] dw, 

JS Jul 

(9) 

where the direction of the surface dS is from inside to 
outside, and the surface integral has been converted to 
a volume integral by Gauss' theorem. The body force 
acting on the fluid within w due to external source is 



/ X(r;t) 

J LO 



duj, 



(10) 



where X is the force per unit volume. The surface force 
acting on the fluid within w due to neighbor fluid is 



7T (r; t) ■ dS 



I 

J LO 



7T diV. 



(11) 



Then as in all hydrodynamical literatures, we have 

R c + R/ = F e +F s , (12) 
which further results in its differential equivalence of 
d 



dt 



[pu] + V r • [puu] = X (r; t) + V r - tt , (13) 



since u> is arbitrary. The last two equivalent equations 
are usually called hydrodynamical equation of motion. 
Actually they are Newton's Second Law for fluids. Later 
we will read them in some views and thus get various 
definitions of the stress accordingly. But before this, we 
have to separate macroscopic motion from microscopic 
motion. 

It is known that each "smallest" element of any fluid 
or material is a particle and the true motion of it is mi- 
croscopic motion. Let us use v (r; t) to denote the mi- 
croscopic velocity of the element. However we can see 
their macroscopic motion only. One example is a flowing 
river, in which we can see the water is transporting, and 
the transporting velocity at each macroscopic point can 



be measured easily, while the motion of each individual 
water molecule is not visible and its microscopic velocity 
can not be easily measured. In most cases, the macro- 
scopic and microscopic velocities are not the same. We 
can also say the macroscopic motion is mass transport, 
and the macroscopic velocity u (p; t) at a macroscopic 
point p is the velocity of mass of centre of the portion 
of the fluid/material in a local small space region Aw p 
around the point. Specifically we have 



u(p;i) / p{r;t)dw 

J Auj„ 



p(r;t)v(r;i)dw, (14) 



where we use the integration representing the correspond- 
ing summation over particles. In this sense, a group of 
microscopic particles can still be regarded as a contin- 
uous medium in a small region. Let us define internal 
velocity as 



w(r;i) = v(r;i) - u(p;i) . 



Then Eq. (14) leads to 



/ 



p (r; t) w (r; t) du; = 0. 



(15) 



(16) 



Let us call the region Aw p as cell p and assume that the 
whole fluid is made of these full cells without overlapping. 
With Eq. (16), the total kinetic energy of the cell p , a 
microscopic kinetic energy, can be written as 



K Aup = \ I p(r;t) v 2 (r; t) dw 

1 JAu n 



.LO D ,mac 



+ K Au 



(17) 



where the corresponding macroscopic kinetic energy 

K A u, p ,mac = t:U 2 (p; t) [ p(r;t)dw, (18) 



and the corresponding internal kinetic energy 



p (r; t) w 2 (r; t) dw. 



(19) 



It is well known that temperature represents internal ki- 
netic energy. As a matter of fact, both macroscopic and 
microscopic velocities and kinetic energies are dependent 
on the velocity of the reference coordinate system. If 
temperature represented microscopic kinetic energy, it 
would also be dependent on the reference coordinate sys- 
tem. On the contrary, the internal velocity and internal 
kinetic energy defined above are independent on the ref- 
erence coordinate system. So is the temperature repre- 
senting the internal kinetic energy. We can not accept a 
temperature dependent on the reference coordinate sys- 
tem. Internal stress should also be independent on the 
reference coordinate system. 

Now let us read every term of Eq. (12), but, for sim- 
plicity, considering the macroscopic motion only. The 
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first left term R c is the total rate of change of momen- 
tum of the system, including the part due to fluid moving 
across the boundary, which is not regular-force-induced. 
Then the second left term R/ is used to cancel this part. 
So the total left side is the total rate of change of regular- 
force-induced momentum. Then the right side must be 
pure regular forces, no imagined force. The external force 
F e is clear and free from problems. The surface force F s 

must be regular forces, so the stress 7r should not have 
the kinetic energy part Tik , and should have the pure 
regular force part 717 only. 

If the term is moved from the left side to the right 
and put into the term F s , the stress becomes 



7T = -pUU+ TTf , 

and Eq. (12) changes into 



R c — F e + 



/ 



dto. 



(20) 



(21) 



For stress definition of Eq. (20), the term puu repre- 
sents the imagined forces between space and space corre- 
sponding to fluid's moving across the boundary, Newton's 
Second Law must be used in the form of Eq. (21), and 
the term R c must include the part due to fluid moving 
across the boundary. So the definition of stress depends 
on how Newton's Second Law being used. 

Next let us read Eq. (12) in view of microscopic mo- 
tion, the complete motion. Then, Eq. (12) becomes 



R c + R/ — F e + F s 



(22) 



with R^ = f u §- t [pv] dw and R^ = f u V r ■ [pvv] du. The 
first left term R^, is the total rate of change of momen- 
tum including the part due to fluid's microscopic moving 
across the boundary, which is not regular-force-induced 
and cancelled by the second left term Rj. Then the right 
side must be pure regular forces according to Newton's 
Second Law, and the stress ir should not have the ki- 
netic energy part. If we move R^ from the left side to 
the right and put it into the stress, then the stress can 
be defined as 



-/9VV+ TTf 



(23) 



with the term — pvv representing the imagined forces be- 
tween space and space corresponding to fluid's micro- 
moving across the boundary, and Eq. (22) changes into 



r; = f 6 + 



/ 



dev. 



(24) 



Let us now separate the microscopic motion into 
macroscopic and internal motions, supposing uj contain- 
ing complete cells only. Then we have 



= / JtIH^ + J; J2 [ P(r;t)w(r;t)duj 
- [pu] dio = R c . (25) 



Actually the rate of flow of momentum is a statistical 
idea as early as in Eq. (9). The first right statistical 
term — pvv of Eq. (23) for a microscopic particle implies 
that the particle has the same probability appearing at 
any point of a certain space region. Let us assume the 
region is the cell in which it resides, then Eq. (16) results 
in the averaged internal momentum at each point 



p(r;t)w(r;t)=0, 



(26) 



which docs not mean no internal motion. Recognizing 
the average nature in the rate of microscopic flow of mo- 
mentum R'j and bringing Eq. (26) into it, we arrive at 



R' f = j> p (r; t) v (r; t) v (r; t) ■ dS 

= £ p (r; t) [u (r; t) + w (r; t)} [u (r; t) + w (r; t)] ■ dS 
= R f + R w , (27) 
where 

R w = <j> p (r; t) w (r; w (r; t) ■ dS. (28) 

Then Newton's Second Law of Eq. (22) becomes 

Re + Rj+R™ =F e + F s . (29) 

As in the above, one has the choice to move R/ from left 

to right and define the tress as n' of Eq. (20) or to keep 
R/ in the left. For either case, we still have a further 
choice whether to put R^ into F s . If we do so but keep 
R/ in the left, the stress becomes Irving and Kirkwood's 
stress [3] 



-p\VW+ TTf 



(30) 



This stress will be zero in the example of figure 4 in 
Zhou's paper [1], as there is no internal motion in the 
material. More interesting is that H w is completely new, 
when Eq. (29) is compared with Eq. (12). The right 
sides of the two equations may be different when mi- 
croscopic motion status is considered in Eq. (29). But 
this difference is only an accuracy problem if it exists. 
The term H w arises from different models of the same 
fluid/material. 

Suppose we study macroscopic properties of a fluid or 
material. As the same macroscopic motion may come 
with different internal motions, the macroscopic proper- 
ties should be averaged ones over internal motions. For a 



3 



given cell with a fixed macroscopic velocity u, if the direc- 
tion of the internal velocity of every microscopic particle 
is rotated with a fixed angle, the macroscopic motion is 
not changed. So let us further assume that the internal 
velocity w for any given microscopic particle has equal 
probability in all directions with the same absolute value. 
Then for any given surface S, we have averaged tensor 
acting on it 

ww-S=-|w| 2 S, (31) 

the averaged internal kinetic energy density 

1 2 3 

£ = 2^l w l =2^ ww ' ( 32 ) 

and finally 

pww = r)kT, (33) 

where the relation e — ^r/kT is used with temperature 
T, Boltzmann constant k, and particle number density 
r\. Please note that pwxv ■ S =f]kTS, the total rate of 
momentum transferred across the surface S, is the net 
rate of that from inside to outside ^r)kT\S and that from 
outside to inside ^rjkT 2 S, so T — \ {T x + T 2 ) even if 
inside temperature T\ and outside temperature T 2 are in 
different cells and not equal. 

Considering this statistical nature, Eq. (28) can be 
simplified as 

R^, = I r/kT (r; t)dS= / fcV p {r)T) doj, (34) 

JS Jul 

which means the net rate of momentum transferred into u) 
due to microscopic particles' internal motion penetrating 
the boundary. At this moment it is an imagined force, 
which can be put inside or outside of stress definition. 
Considering the internal motion condition of Eq. (16), 
we can easily recognize that the increased momentum 
due to microscopic particles' internal moving across the 
boundary should contribute to the macroscopic motion 
in lu, described with u (p; t). 

Now let us consider measurements of the fluid. 
When some detector is placed into the fluid to de- 
tect/measure some property, the microscopic particles 
will inevitably collide with it and the collision must be 
detected/measured. If the detector is not moving with 
the same velocity of the fluid's macroscopic motion, the 
collision will be complicated and depend on many fac- 
tors. So a point detector moving with the same macro- 
scopic velocity of the fluid at the point can measure some 
intrinsic-like properties. In this case, the collisions be- 
tween the fluid particles and the detector is due to par- 
ticle's internal motion only, and the averaged force in 
the collisions can be easily proven as pww • AS —rjkT AS 
for the detector's small surface area AS exposed to the 
particles, supposing the collisions are complete elastic. 
So we found some physical force background for R^, or 



fcVp (r/T). If we regard it as an additional regular force, 
we can define the internal stress as 

^ mt =-kf]T I +^c~f, (35) 

with unit tensor I , while Newton's Second Law be- 
comes 

R c + R/ = F e + Fj„ t , (36) 

where 

Fj nt = f 7r int (r:t) -dS = / V r - 7T mt dto. (37) 

JS Jul 

Essentially n int is the same as the stress n in our 
previous work [2] , where the lattice translation motion is 
macroscopic and the motion of particle in the MD cell is 
internal for every cell. 

As temperature is a fundamental physics quantity, we 
may also write Newton's Second Law as 

R c + R / = F e + F s + F t , (38) 

where the temperature force F t = — R w . A simple ex- 
periment can be done as follows to see the action of this 
force. Just have a kettle with half-filled water, boil it a 
few times to get rid of bubbles, cool it and let it calm 
down to no visible motion for a while, and heat it again. 
Then the water flow is visible even without bubble and 
before boiling. This flow is due to the temperature force, 
as the temperature is not even. 

In summary, we have discussed many possibilities of 
stress definitions for different models of the fluid/material 
and for different forms of using Newton's Second Law. 
However any stress definition will be meaningful only 
when it is used in Newton's Second Law and Newton's 
Second Law is used correctly. By separating macroscopic 
motion and internal motion from microscopic motion, we 
found that the temperature force can also drive a fluid 
flow. 
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